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Euclidean space hadronic correlators are computed in quenched QCD at small quark mass using 
truncations of quark propagators which include or exclude low eigenvalue eigenmodes of the Dirac 
operator. High modes provide the dominant contribution to parity averaged correlators, especially 
at short distances. Differences of correlators of opposite parity receive most of their contributions 
from low modes and are much smaller in size than parity averages at short distances. The pion 
propagator in any correlator to which it couples receives a large contribution from low modes, while 
the tensor meson correlator receives a tiny contribution from low eigenmodes. 

I. INTRODUCTION 

Are all hadrons alike? This question was asked most famously jjj in the context of QCD sum rules, but it appears 
repeatedly in phcnomenological descriptions of hadronic structure. In this work I investigate a particular kind of 
similarity: how eigenmodes of the 4-dimcnsional Dirac operator contribute to Euclidean space hadronic correlators, 
and so (indirectly) how these eigenmodes contribute to hadronic spectroscopy and couplings in different channels. My 
results might have implications for the questions of whether excited states of mesons show features of chiral symmetry 
restoration (namely, that the masses of various states, and their couplings to currents, show degeneracies), what the 
origin of these features might be, and how excited states differ from low lying meson states. 

Lattice simulations encode "snapshots" of the QCD vacuum, and could presumably address questions about the 
composition of hadronic states. This subject has been investigated before @,G|. The idea is to compare the behavior of 
hadronic correlators built from quark propagators which correspond to a truncated set of eigenmodes of the massless 
Dirac operator with that of correlators computed "exactly." Low eigenmodes of the Dirac operator saturate the 
pseudoscalar and axial vector correlators at large distance and contribute somewhat to the nucleon channel, less so 
the vector and Delta channels. They also make the dominant contribution to the eta-prime channel, as expected 
from the Witten-Veneziano formula jjj]. In this work I revisit these studies, with an eye toward contrasting the 
Dirac eigenmode composition of the long-distance part of hadronic correlators with the short-distance parts, as well 
as the composition of correlators which are sums and differences of parity partners (75 and 1, 7^ and 7475). All 
these calculations reveal the following qualitative features of mesons in quenched QCD: Low eigenmodes of the Dirac 
operator do not affect the part of correlators where high-lying states appear. These low modes are the ones which 
determine the quark condensate via the Banks-Casher [5j relation, as well as eigenmodes "at the QCD scale" (a few 
hundred MeV) . Other ways of describing the properties of low modes are that they are strongly influenced by chirality- 
mixing interactions, and that they are the main contributors to the strong interaction observed in the pseudoscalar 
current-current correlator. 

High modes contribute strongly to hadron correlators at short distance, where the correlators would be dominated 
by excited states. They contribute less to the long distance part of meson correlators. They do not make a major 
contribution to parity asymmetric correlators. In this sense, valence quarks in high-lying states decouple from the 
condensate. This also means that even for light quarks, chiral symmetry breaking and confinement are not directly 
related. These observations are rather simple, and perhaps one does not "need" lattice simulations to see them. I 
suspect that no phenomenologist will find them surprising (although I suspect that different phenomenologists will 
be unsurprised for different reasons.) Instanton models plus lattice simulations have also been used to compute the 
amplitude for chiral mixing of qq correlators @, Q • Ref . also observed that the lowest eigenvalue modes had a 
large autocorrelation in both density and chirality at short distance, which slowly dies away as the eigenmode rises. 

Analyses of hadronic spectroscopy H, ||J 0, Ujl 0, 0, Q have been used to argue that excited states of the 
meson and baryon spectrum are effectively chirally restored, in the sense that particles of opposite parity are nearly 
degenerate in mass. It is argued that this "effective restoration of chiral symmetry" is a natural consequence of the 
soft nature of chiral symmetry breaking: generally in quantum mechanics, low excitations are sensitive to symmetry 
breaking but high excitations are not. For mesons, the crossover to chiral symmetry is argued to occur in the range 
1.5-2 GeV. 

In perturbation theory, correlators of currents which are chiral partners of each other (the vector and axial vector 
currents, for example) are identical, because QCD is a vector theory and chirality is conserved at all quark-gluon 
vertices. Condensates allow one to parameterize nonperturbative physics, and can lead to differences in correlators. 
However, QCD is a confining theory, and the spectrum of QCD is one of bound states, so identity of the correlators 
must have consequences for spectroscopy and matrix elements. That is especially true for the large- N c limit of QCD, 
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where all excitations are narrow resonances. Arguments based on semi-local duality favoring effective chiral restoration 
at high excitation have been made but have been criticised as being cutoff dependent 0]. 

A summary of ideas supporting chiral restoration and its connection to the quark model has recently been given 
by Swanson|17|. 

In instanton models |l8j low eigenmodes of the Dirac operator are built from an overlap of single- instanton zero 
modes, and are strongly chirally asymmetric. As the magnitude of the quark eigenmode rises, the eigenfunctions 
couple less and less to instantons, and there should be a "crossover" to chirally-symmetric physics. 

The common lattice correlators used to extract spectroscopy involve projections onto momentum eigenstates. For 
example, the p = correlator is found by by averaging over time slices: 

Cj{t) = 52(0\Oj(x, «)O i (0,0)|0) 

n 

(i) 

A second kind of observable is just the point-to-point correlator 

Ui(x) = Tr(J?(x) J?(0)), (2) 

where the current will be proportional to Jf{x) — xjj(x)T a T(i)ip(x), and its ratio to the free-field correlator Ri(x) — 
Hi(x) /U®(x). It is very difficult to compute masses of highly excited states from lattice simulations. Their signals 
vanish exponentially compared to (and underneath) the lightest state in the channel. Each new state requires a fit 
with two more parameters (the mass and (0|Oj|n)| 2 ), which in turn requires a fine lattice with many lattice points. 
At very short distance (order one lattice spacing) lattice discretization artifacts affect results. (Note that correlators 
like C(t) probe long distances, even at small t.) I have not done direct calculations of spectra, and so my approach 
is more indirect and my results are tentative. 

I will look at correlators with point sources and sinks, Oj(x,t) = ip(x,t)Tjtp(x,t). Lattice simulations usually do 
not use pointlike currents as interpolating fields because they do not couple to low states as well as more extended 
operators, and parameters of the lowest states are usually the goals of the simulation. However, for our purposes, it 
will be interesting to keep this simple form, and look at I^ 's which are chiral partners. 

In contrast to usual lattice simulations, all results presented here are qualitative. 



II. EXACT RESULTS FROM OVERLAP ACTIONS 

The calculations presented here are done with overlap [l^ fermions. This lattice discretization preserves a lattice 
version of exact chiral symmetry at nonzero lattice spacing, without flavor doubling, making it particularly useful for 
addressing questions associated with chiral symmetry. The following properties of overlap actions are relevant to this 
work: 

The eigenmodes of any massless overlap operator are located on a circle in the complex plane of radius xq with a 
center at the point (xo, 0). The corresponding eigenfunctions are either chiral (for the eigenmodes with real eigenvalues 
located at A = or A = 2xq) or nonchiral and paired; the two eigenvalues of the paired nonchiral modes are complex 
conjugates. 

The massive overlap Dirac operator for bare quark mass m is conventionally defined to be 

Tfl 

D(m) = (l-—)D(0)+m (3) 
2x 

and it is also conventional to define the propagator so that the chiral modes at A = 2xq are projected out, 

D- 1 (m) = l — -m-\m)-—). (4) 

v ' l-m/{2x ) y y ' 2x ' 

The contribution to the propagator of a single (positive chirality) zero mode in the basis where 75 = diag(l, —1) is 

i>e»)-' = = ( l3 ' +> 0+l °V o» 



in 
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Nonzero eigenvalue eigenmodes of D(0)t_D(0) are also chirality eigenmodes. The jth pair of nonchiral modes con- 
tributes a term to the propagator 



where, using D(0) 1 'D(0)|j7i} = \j\jh) for chirality h, fi = m,/(2x ), and ej = Xj/(2xo), the entries are 



v i-s 2 ) 

- 2^ £ 2 +>(l- e 2 )' (8) 



(The eigenmodes of -D(O) have eigenvalues 2xo(e 2 ±iej^J 1 — e 2 ).) For a summary of these (and other) useful formulas, 

see Ref. (2^ (for the special case xo = 1/2). 

What is important about these well-known results is their application to correlators expressing the sum and differ- 
ence of opposite parity channels: 

C± r (x,x') oc Tr£>- 1 (m)rA±r (9) 

where 

A± =i)- 1 (m)T75^ 1 M75 (10) 

The parity difference is particularly simple: it involves the sum combination, which because of the Ginsparg- Wilson 
relation [21j is equal to 

A_ = 2mlD- 1 (m)^D- 1 (m) (11) 

The volume-summed trace of A_ gives the Gell-Mann-Oakes-Renner (GMOR) relation. Eigenmode by eigenmode, 
A + and A_ can be evaluated directly from Eq. the former contribution comes from the off-diagonal part of D^ 1 
while the latter contribution comes from the diagonal elements. 

This decomposition makes no reference to any dynamics. Thus the parameter which differentiates between chiral 
symmetry and asymmetry is just the quark mass-it is the only chiral breaking parameter which is available. The ratio 
of a contribution of an eigenmode to the chiral-difference correlator to the chiral symmetric correlator scales roughly 
as m/Xj. This is rather similar to continuum free-field behavior, with A_(p) = m/p 2 , A + (p) = p/p 2 . It is NOT the 
behavior seen for Wilson- like fermions, which have an explicit chiral-symmetry breaking term at large momentum, 
D~ 1 (m) = m + i ^2 In sin(p A1 a)/a + 2a sin (p^a/2). Staggered fermion correlators entangle parity partners in a 
way I do not know now to present simply. 

Naively, this means that, unless the to — > limit is singular, parity difference correlators will go to zero with the 
quark mass. However, we have to be more careful j23j, because the eigenmode density p(\) diverges in the ultraviolet 
as A 3 . Mode sums generally require subtraction: high eigenmodes do contribute to all susceptibilities. For example, 
the pseudoscalar susceptibility (time integral of Cj(t) for Tj =75) is 

X" = V E^K(f)> - / dX P^y^ (12) 

x,y 

ao = \ $>S (*K(y)> - / dX P (X) A 2 2 " ™' (13) 



and the scalar susceptibility is 



x.y 



For either susceptibility, it is necessary to break the integral over A into two parts, one for the low modes < 
A < A c where A c >> to, and one for the high modes. The low modes give the Banks-Casher relation for 
Xtv = p(0)/(nm) + . . .. The high frequency part of the integral must be subtracted twice, 



1 " A 



TO 



X*=—(l dX P( X )jTT — 2) +7o + 72m 2 +TO 4 ^(A c ,77i). (14) 
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The scalar susceptibility has a similar behavior, except that there is no contribution from A = |22(. For overlap 
fcrmions, the GMOR relation is exact, 



2 - 

X* = —W${m)) (15) 



m 



and the scalar susceptibility is 



Xao =2^-(^(m)), (16) 
am 



so that the expansion coefficients in the expansions for x-rr an d Xa are related. In the difference x-k — Xa the 
contribution from the chiral condensate (which comes from eigenmodes near A = 0) survives and dominates in the 
772 — ► limit. The mass-independent term in Xn — Xa cancels. The chiral sum Xn + X.a retains the high-A dependent 
7o coefficient. The terms proportional to positive powers of m vanish in the massless limit in either case. What is 
true for an integral is not necessarily true for an integrand, but it is a plausible assumption to expect that chiral 
difference correlators themselves will receive large contributions only from low eigenmodes, while high eigenmodes 
will contribute strongly to chiral sum correlators. For example, one might expect that the low modes would dominate 
the pseudoscalar correlator: then 0, 

C PS (x,x') = ({ A, x\j, A, x) * {j, A, x'\j, A, x')}) (17) 

\ Z + 777 Z 

This reasoning is dangerous because of the UV divergence of p(A), but it can be tested by simulation. 

Zero modes decouple from C+(x,x'). They make a large contribution to C-(x,x'). This is not surprising: zero 
modes are chiral. They will, however, considerably distort chiral-difference correlators in small volumes. 



III. EXAMPLES 

I have done a set of quenched spectroscopy runs using a particular implementation of the overlap operator j2^| 
with HYP (2(| gauge connections. Eigenmodes of the Dirac operator are computed using the conjugate gradient 
algorithm of Ref. [2jJ. These studies are done on a small data set of quenched configurations. It has 20 16 4 lattices 
generated with the Wilson gauge action at coupling (i — 6.1, corresponding to a lattice spacing of about 0.09 fm. 
These lattices have the same lattice spacing and quark masses as were used in a large scale matrix element simulation 
pEjj. so all necessary hadron masses and matrix elements are in principle known. I have computed the lowest 20 
eigenmodes of Z?(0)^I?(0) and recoupled them into eigenmodes of D: there are N chiral zero modes and 2(20 — iV ) 
nonchiral paired modes. Eigenmodes of the overlap Dirac operator whose eigenvalue is less than around 500 MeV 
have been computed. The spectrum of the imaginary part of the eigenvalue is shown in Fig. ^ The energy is inferred 
from the lattice data assuming an inverse lattice spacing of 2200 MeV. I compute hadron correlators using full quark 
propagators, and several kinds of quark propagators built of truncated mode sums: quark propagators with only zero 
modes, propagators without zero modes, and propagators from which the lowest 10 or 20 eigenmodes of D^D have 
been excluded. 

I will show results corresponding to bare quark masses in lattice units of am q = 0.020, which is about half the 
strange quark mass. A meson made of a pair of am q = 0.020 quarks is an approximation to the physical kaon. Results 
for other small quark masses are similar. At higher quark masses (above am q = 0.050 or ttt^ > 2.5m s ) the qualitative 
features I present do no longer occur. Pictures shown in Refs. 0, Q show that low modes cease to contribute a 
dominant part of any correlator. 

Because these simulations are done in finite volume, and in quenched approximation, some channels have contribu- 
tions from zero modes of the Dirac operator. These contributions would not be present in infinite volume. In other 
channels these modes arc absent. 



A. "Spectroscopic" correlators with point sources 



I begin first with correlators C (t) of Eq. ^ Results for pseudoscalars and scalars are shown in Fig. [3 Panel (a) 
shows the sum of pseudoscalar and scalar correlators. Note the complete saturation of the correlator for t > 5 by 
low modes, and the essentially complete saturation of the correlator by high modes at shorter distances. The pionic 
contribution to the correlator, using (0|t/>75 ip\PS) = mp s j p 5 / '(2m) , (with the fitted parameters mps and fps from 
Ref. 28]) is superimposed. The low modes make up the bulk of the pion's contribution to the correlator. 
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FIG. 1: Histogram of eigenvalues of the Dirac operator of eigenmodes recorded for this study. The peak at the lowest box 
contains all zero modes. The black solid lines curve shows the distribution of the twenty lowest eigenmodes of D< D, while the 
red dashed lines separate the distribution of ten modes. 

Panel (b) shows the difference of pseudoscalar and scalar correlators. This signal is heavily influenced by zero 
modes, so that while low modes again contribute a "pion-like" signal, most of this contribution involves at least one 
zero mode. High eigenmodes make a negligible contribution to the correlator at any time separation. 

The axial-scalar (7075 — 7570) ± (70 — 7o) combination shows similar behavior (panels (c) and (d). The axial vector 
current decouples from the Goldstone mode in the chiral limit, and this is reflected in the "sum" correlator by a slower 
saturation of the pion contribution to the correlator by low modes. Again, the bulk of the signal at small time steps 
comes from high modes. High modes do not contribute to the chiral difference. 

Vector and axial combinations (Fig. |3J panels (a) and (b)) lack any single state which is dominated by low modes. 
High modes again make a tiny contribution to the chiral-asymmetric channel, and saturate the signal in the chirally 
symmetric channel. The absolute size of the correlators in the two channels is quite different at short distances. 

The correlator C+(t) is almost two orders of magnitude larger than C-(t) in the vector-axial vector channel, and 
almost three orders of magnitude larger in the pseudoscalar-scalar channel. This is not a direct observation that the 
spectrum of quenched QCD is parity-doubled at high excitation. However, this is the "raw data" from which masses 
are extracted, so it is likely that any fit to it for masses will produce parity symmetric results. 

Perhaps direct comparisons, in Fig 01 are clearer. Panel (a) shows the pseudoscalar and scalar correlators; I have 
subtracted the measured pion contribution from the pseudoscalar channel. Panel (b) shows the vector and axial vector 
correlators. The correlators Cj(t) defined in Eq. ^ where Tj are parity partners (75 and 1, 7$ and 7^75), become 
equal for t less than about 0.2 fm. 
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FIG. 2: Combinations of point-to-point pseudoscalar and scalar (n ± ao) sum (a) and difference (b) correlators, and axial 
vector 7075 — 7075 and scalar 70 — 70 sums (c) and differences(d). The full correlators are shown by open squares. Thheir 
approximation by restricted-mode quark propagators with 20 lowest modes (red octagons) or 10 lowest modes (green bursts). 
Crosses and diamonds show the contribution to the correlator from quark propagators with 20 and 10 low modes excluded. 
The curve is the contribution to the correlator from the pion, using its measured mass and decay constant from a different 
simulation. In panels (b) and (d) the contribution of the lowest 20 eigenmodes with zero modes excluded is shown in blue fancy 
diamonds, and the pure zero mode contribution (magenta fancy crosses) is shown in (b). 



B. Point-to-point correlators 

One can also compute point-to-point correlators (Eq. The short distance, parity-summed correlators receive 
little contribution from low eigenmodes. That their ratios approach unity at small x is just asymptotic freedom. 
However, this is a confining theory: the spectrum of quenched QCD consists of zero width resonances plus quenched 
artifacts (mostly associated with eta prime hairpins). A more correct statement of what is seen is that the resonances 
which make up the low-x part of the sum correlators receive little contribution from low modes. Results are shown in 
Fig. El 

The parity difference correlators vanish at low x. Low eigenmodes (including zero modes) dominate the correlators 
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FIG. 3: Point-to-point (7, — 7*) ± (7175 — 7i7s) (rho plus/minus ai) correlator: (a) sum, (b) difference, and (c) tensor or 
a,2 correlator (7*7? — 7»7j). Full correlators are shown as black squares, and their approximation by restricted-mode quark 
propagators with 20 lowest modes (red octagons), fO lowest modes (green bursts), and pure zero mode contribution (magenta 
fancy crosses) are also shown. The contribution of the lowest 20 eigenmodes with zero modes excluded are blue fancy diamonds. 
Crosses and diamonds show the contribution to the correlator from quark propagators with 20 and 10 low modes excluded. 



at large x, where the lightest states in the channel appear. Alternatively, the strong attraction seen at large x in the 
pseudoscalar channel comes almost entirely from eigenmodes below 500 MeV. 

In the absence of reliable spectroscopy calculations, we can roughly quantify the "breakpoint" between low and high 
mode contributions by considering the contribution of a single resonance of mass m to one of these ratios. This curve 
is proportional to the ratio x 5 K(m, x), where K(m, x) is the free field propagator for a particle of mass m. A family of 
curves of varying masses (in lattice units) is shown in Fig. [S] the height has been rescaled as m 2 to produce a plateau 
of maxima. What is important from this picture is not the height of the peaks, it is their location. If we compare 
the high eigenmode and full propagator curves in FigOJa) and (c), we see that they separate at a lattice distance 
of about 5 units. Any resonance with a lattice mass lighter than amn < 0.75 would make its peak contribution at 
larger x. This mass is about 1700 MeV, and is the rough dividing point between hadrons built dominantly of chiral 
sensitive modes and ones which are not, for quenched QCD. (Again, recall that the inverse lattice spacing of these 
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FIG. 4: (a) Point-to-point (75 — 75) correlator, with pionic contribution removed (squares) and 1 — 1 (crosses) correlators at 
bare mass am q = 0.02. (b) point-to-point (7; — 7i) (squares) and (7*75 — 7*75) (crosses) correlators 

simulations, above which discretization effects dominate physics, is about 2.2 GeV). Swanson's estimate H3 of 2.5 
GeV for the crossover corresponds to ma ~ 1.1, but given the roughness of either estimate, I would not take the 
difference seriously. 

Finally, we show in Fig. [7]thc saturation of the amplitude for chirality flipping proposed by Ref. [2j| and computed 
on the lattice in Ref. la, 



U„(x)+Il ao (x)' 



(18) 



using complete quark propagators and high-eigenmode truncations. Clearly the high eigenmodes make a tiny contri- 
bution to R , it is the low modes that are strongly influenced by chirality mixing interactions. 



IV. CONCLUSIONS 



Lattice simulations in which contributions to quark propagators from low and high eigenvalue eigenmodes of the 
Dirac operator are separated identify the following qualitative features of light-quark mesons in quenched QCD: Low 
eigenmodes make a large contribution to the pion propagator. They are responsible for the strong long-distance at- 
tractive interaction seen in the pseudoscalar channel. Low eigenmodes of the Dirac operator make a small contribution 
to the short distance part of correlators (where excited states contribute). These are the eigenmodes which determine 
the quark condensate via the Banks-Casher relation, as well as eigenmodes at the QCD scale" (a few hundred MeV). 
This implies that valence quarks in high-lying states decouple from the condensate. 

I have not directly observed parity doubling at high excitation in the meson spectrum, but lattice correlation 
functions in parity-partner channels (from which fits to masses would be performed) are equal at a part in 10 3-4 at 
short distances. Point-to-point correlators suggest that hadrons of mass above 1.7 GeV are insensitive to low Dirac 
eigenmodes. Low modes do not contribute to the tensor meson correlator. 

Perhaps these results can be used to constrain phcnomenological models of hadron structure. Results such as Fig. 
EJa) suggest that it might be profitable for lattice calculations done at small quark mass to exploit low eigenmodes of 
the Dirac operator in simulations. 
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FIG. 5: Point-to-point (7 M — 7^) ± (7 M 7s — 7 M 75) vector-axial vector sum (a) and difference (b), and (75 — 75) ± (1 — 1) 
pseudoscalar-scalar sum (c) and difference (d), normalized by free field vector and scalar currents, Full correlators are shown 
in black and approximation by restricted-mode quark propagators with 20 lowest modes (red crosses) and contributions from 
propagators containing all but the lowest 20 (blue diamonds) or 10 (green squares) modes. 
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